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^Jj Abstract 

The energy conditions are derived in the context of /(-R, T) gravity, 
where R is the Ricci scalar and T is the trace of the energy-momentum 
tensor, which can reduce to the well-known conditions in f{R) grav- 
ity and general relativity. We present the general inequalities set by 
the energy conditions in terms of Hubble, deceleration, jerk and snap 
parameters. In this study, we concentrate on two particular models 
of f{R,T) gravity namely, f{R) + XT and R + 2f{T). The exact 
power-law solutions are obtained for these two cases in homogeneous 
and isotropic f{R,T) cosmology. Finally, we find certain constraints 
which have to be satisfied to ensure that power law solutions may be 
stable and match the bounds prescribed by the energy conditions. 
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1 Introduction 

Recent astrophysical observations form supernova type la}\ cosmic microwave 
background anisotropies^-* , large scale structure'^^ baryon acoustic oscillations^ 
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and weak lensing^^ indicate that the universe is accelerating in the current 
epoch. The most promising feature of the imiverse is the dominance of exotic 
energy component with large negative pressure, known as dark energy (DE). 
A number of alternative models have been proposed in the framework of gen- 
eral relativity (GR) to explain the role of DE in the present cosmic accelera- 
tion. Unfortunately, up to now, no suitable candidate is found, which boosts 
our interest in modified theories of gravity. Firstly, the Einstein-Hilbert 
action has been modified by replacing scalar curvature R by an arbitrary 
function of i?, this theory is known as f{R) gravity^-*. The other alternative 
theories of gravity include /(T) gravity^^ where "7^' is the torsion scalar in 
teleparallel gravity, Gauss-Bonnet gravity^) and f{R,T) gravity^^ 

The f{R,T) gravity is the generalization of f(R) gravity involving the 
dependence of the trace of energy-momentum tensor T. The dependence 
of T may be induced by exotic imperfect fluids or quantum effects. The 
cosmological reconstruction of f{R,T) gravity has been studied in recent 
literature^" In a paper^\ the reconstruction of FRW cosmology is pre- 
sented for f{R,T) = R + 2f(T) model. Houndjo and Piattella^^-* constructed 
f{R, T) models describing the uniflcation as well as transition of matter dom- 
inated phase to late accelerating phase. The chaplygin gas f{R,T) models 
are investigated in^^~^^) and it is shown that dust fluid reproduces ACDM, 
Einstein static universe and phantom cosmology^^^ In our previous work^^\ 
we have reconstructed some explicit models of f{R, T) gravity for anisotropic 
universe and explored the phantom era of dark energy. We have also discussed 
the validity of flrst and second laws of thermodynamics in this modifled 
gravity^'^^ The existence of exact power law solutions for FRW spacetime 
has been investigated in modifled theories of gravity Here, we shall 
show that FRW power law solutions exist for a particular class of f{R,T) 
gravity. 

The classical energy conditions of GR are profound to the Hawking- 
Penrose singularity theorems and classical black hole laws of thermodynamics^^) . 
These conditions have been used to address several important issues in GR 
and cosmology^*^-'. Many authors have investigated the energy conditions in 
the context of modified theories including f{R) gravity^^^^*^\ f{R) gravity 
with nonminimal coupling to matter^^-*, modified Gauss-Bonnet gravity^^\ 
modifled f{G) gravity with curvature- matter coupling^^\ Brans-Dicke theory^^\ 
/(T) gravity^^) and f{R,T) gravity^^^ In recent work^^\ the authors study 
the energy conditions for a special form of f{R,T) gravity, f{R,T) — R + 
2/(r) and discussed the stability of two f{T) models. 
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In this work, we are interested to set energy conditions bounds on exact 
power law solutions in f{R, T) gravity. The FRW power law solutions are ob- 



tained for f{R, T) = f{R) + AT and f{R, T) = R + 2f{T) gravity. We derive 



the energy conditions for more general as well as particular class of f{R,T) 
gravity. The standard form of energy conditions in GR and f{R) gravity 
can be recovered in the limit of f{R,T) — f{R) and f(R,T) = R. We show 
that for f{R,T) = f{R) + AT, the null energy condition (NEC) and strong 
energy condition (SEC) can be derived by using the Raychaudhuri equation 
with the requirement that gravity is attractive. The resulting inequalities for 
NEC and SEC are equivalent to the energy conditions obtained in terms of 
effective energy-momentum tensor. 

The paper is organized as follows: In next section, we present the general 
formulation of the field equations of f{R,T) gravity in FRW cosmology. In 
section 3, the energy conditions are derived and hence presented in terms of 
deceleration {q), jerk (j) and snap (s) parameters. Section 4 is devoted to 
obtain exact power law solutions for two specific forms of /(i?, T) gravity. 
We also analyze the constraints of energy conditions for these models. In 
section 5, we investigate the perturbation and stability of power law solutions. 
Finally, section 6 summarizes the obtained results. 



The f{R, T) theory of gravity is an interesting modification to the Einstein 
gravity by introducing an arbitrary function of scalar curvature R and trace 
of the energy-momentum tensor T. The action for this theory coupled with 
matter Lagrangian C{rnatter) is given by^^ 



where g is the determinant of the metric tensor g^^^, we use the units SttG = 
c—1. The energy- momentum tensor of matter is defined as^^^ 



2 f{R,T) Gravity 




(1) 




(2) 
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Varying this action with respect to the metric tensor, we obtain the field 
equations of f{R,T) gravity as 

= T^u-fT{R,T)T^u-MR,T)Q^,, (3) 

where fn^R.T) and fx^R^T) denote derivatives of f{R,T) with respect to 
R and T respectively; □ = g^'^V ^ is the d'Alembert operator, is the 
covariant derivative associated with the Levi-Civita connection of the metric 
tensor and is defined by 



The contribution to the energy momentum tensor of matter is defined as 

T^,y = {p + p)u^Uy - pg^^, 

where is the four velocity of the fiuid, p and p denote the energy density 
and pressure, respectively. We can take C{matter) = —p, then 9^,^ becomes 

= -2T^^ - pg^^. (5) 

Consequently, the field equations ([3]) can be expressed as effective Einstein 
field equations of the form 

R,u - \r9,u = T^JJ, (6) 

where T^j^-^ is the effective energy-momentum tensor in f{R, T) gravity which 
is defined as 



1 + MR, T))T,, + pg,jT{R, T) + T) 



- RfdR, T))g,, + (V^V. - g,u^)fR{R, T)] . (7) 
We consider the homogeneous and isotropic fiat FRW spacetime as 

ds^ = dt^ - a^{t)dx^, (8) 
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where a{t) is the scale factor and dx^ contains the spatial part of the metric. 
In the FRW background, the field equations may be rewritten as 



(9) 



where pe// and Pejf are the energy density and pressure respectively, defined 
as 

1 



Peff = 

Peff = 



fn 
1 

fn 



P+{P + P)fT + 2 (/ - ^/«) - ^H{RfRR + T/rt) 
p - i(/ - RfR) + {R + 2RH)fRR + R^fRRR 



+ 2RffRRT + {f+2fH)fRT + t''fRTT 



, (10) 



(11) 



the Hubble parameter H is defined hy H = a/a and dot denotes derivative 
with respect to cosmic time t. 



3 Energy Conditions 

Raychaudhuri equation is the key to SEC and NEC together with the re- 
quirement that gravity is attractive for a spacetime manifold endowed with 
a metric g^u^'^^- For the congruence of timelike geodesies defined by vector 
field u^, the Raychaudhuri equation reads, 

(if) 1 

J? = - ^'^'^^/^^ + ^'^'^^/^^ - ^/^-^''^^ (12) 

where is the Ricci tensor, and 9, a^y and uj^y are the expansion param- 
eter, the shear and the rotation associated with the congruence respectively. 
The evolution equation for 9, the expansion scalar of a congruence of null 
geodesies defined by the null vector field k^, is given by 

JQ 1 

= --^2 _ ^/..^^^ ^ _ i^^M/,-. (13) 

or z 

Raychaudhuri equation is known to be purely geometric and hence, develops 

no reference to any theory of gravity. As the shear tensor is purely spatial 
(^^^(^nu ^ 0, thus, for any hypersurface of orthogonal congruence (o;^,/ = 0), 
the conditions for attractive gravity become 

SEC : Rf,yuV ^ 0, NEC : R^^k^k" ^ 0. (14) 
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One can use the field equations of any gravity to relate i?^,^ to the 
energy- momentum tensor T^^,. Thus, the combination of the field equations 
and Raychaudhuri equations can set the physical conditions for the energy- 
momentum tensor. In the framework of GR, the conditions f lT^ can be 
written as 

Rf^un'^u" = (T^, - '^g^,)u''u'' ^ 0, R^^k'^k" = T^.k^k' ^ 0. (15) 

For perfect fiuid, this equation is reduced to the well-known form of SEC and 
NEC in GR, 

p + 3j9^0, p + p^O. (16) 

In modified theories of gravity including f{R) and /(T) gravity, R^y can be 
obtained in terms of the energy-momentum tensor by using the corresponding 
field equations. However, this does not seem apparent in f{R,T) gravity. 

We consider the effective energy-momentum tensor T^^-^, so that the con- 
ditions in Raychaudhuri equations are represented as 

Hence, the energy conditions in GR can be applied by replacing energy den- 
sity p and pressure p by pejf and Peff, respectively. Since the Raychaudhuri 
equation holds for any geometrical theory of gravity, we will keep the phys- 
ical motivation of focussing of geodesic congruences along with attractive 
property of gravity to develop the energy conditions in f{R,T) gravity. We 
also assume that standard matter obey the four energy conditions. Using 
the effective modified field equations ([6]), the energy conditions for f{R,T) 
gravity are given by 

NEC : 

Peff + Peff = 

+ 

WEC : 

Peff = 
+ 



1 

Ir 



(p + p)(l + /t) + (i? - RH)fRn + R^fRRR 



'^RTfRRT + (T — TH)fRT + T'^Irtt 



^ 0, 



(17) 



1 

Ir 
TJrt 



p + {p+p)fT + -if - rir) - mRfRR 



> 0, Peff + Peff > 0, 



(18) 
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SEC : 

Peff + ^Peff 



ip + 3p) + (p + p)fT -f + RfR + ^R'fRRR 



1 

Ir 

+ 3{R + RH)fnR + 6RffRRT + 3{f + fH)fRT 

+ 3f^fRTT\>0, Peff + Peff >0, (19) 



DEC : 

Peff - Peff 



'RRR 



^ \{p-p) + {p + p)fT + f-RfR-R'ff 

JR L 

- {R + 5RH)fRR-2RffRRT-{f + 5fH)fRT 

- f^fRTr] ^ 0, Peff + Peff > 0, Peff ^ 0. (20) 

The inequalities f|T7j) - fl20l) represent the null, weak, strong and dominant 
energy conditions in the context of f{R,T) gravity for FRW spacetime. 

We define the Ricci scalar and its derivatives in terms of deceleration, 
jerk and snap parameters as 

i?= -6/7^(1 -g), R= -6H\j -q-2), 

R= -QH\s + q'' + 8q + Q), (21) 



where^*^ 



Id la 1 "a 

i/^ a H-^ a a 



Using the above definitions, the energy conditions ( IT711 - (l20l) can be rewritten 
as 

NEC : {p + v){l + fT)-QH\s-3 + {q + l){q + 8))fRR + H^mU 

- q - 2)ffRRR - 12H\3 - q - 2)ffRRT + (f - fH)fRT + f 

X fRTT ^ 0, 

WEC : p + {p + p)fT + \f + ?>H\l-q)tR + l9.H\3-q-'^)tRR 

- 3HffRT)^0, Peff + Peff 

SEC : {p + 3p) + {p + p)fT-f-6H\l-q)fR + 3[6H%j-q-2)f 

X fRRR - I8{s +j + q' + 7q + A)fRR - 36H%j - q - 2)ffRRT 

+ 3{f + fH)fRT + 3f^fRTT>0, Peff + Peff 
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DEC : {p-p) + {p + p)fT + f + 6H\l-q)fn-[6H%j-q-2)]^ 
X fRRR - 6H\s + 5j + (g - l)(g + 4))/^^ + 12H%j - q - 2)f 

X fuRT - (T + 5tH)fRT - f^fRTT > 0, Peff + Pef f > 0, 
peff > 0. 

The energy conditions in /(-R) gravity^^-* can be recovered for /(-R, T) = f[K) 
and also in case of GR for particular choice f{R,T) = R. To illustrate how 
above conditions can be used to place bounds on f{R, T) gravity, we consider 
two particular forms of f{R,T) gravity, 
(i) /(i?) + AT, (ii) i? + 2/(T). 

We shall obtain the power law solutions for each case and hence the 
constraints set by the respective energy conditions. 



4 Power Law Solutions for /(i?, T) Gravity 

It is important to study the existence of exact power solutions corresponding 
to different phases of cosmic evolution. Such solutions are particularly rel- 
evant because in FRW background they represent all possible cosmological 
evolutions such as radiation dominated, matter dominated or dark energy 
eras. We discuss power law solutions for two particular models of f{R,T) 
gravity. 

4.1 f{R,T) = f{R) + XT 

For the particular case f{R, T) = f{R) + XT^^'^^\ the effective Einstein field 
equations are given by Eq.([n]) with 



^eff ^ _±_ 
f. 



R 



(22) 



where Jr is the derivative of f{R) with respect to scalar curvature R. The 
Friedmann equation and the trace of the field equations are given by 



02 



3 

Ir 



XT 1 
p+A(p + p) + — + 



&RI 



RR 



2 2 

RfR + 3nfn{R, T) - 2/ = (1 + 3A)r + 4Xp, 



(23) 
(24) 
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where = 3d/a is the expansion scalar. 

The standard matter satisfies the following energy conservation equation 

p = -0(p + p). (25) 

For the homogeneous and isotropic spacetime, the field equations can be 
represented by Raychaudhuri equation 



+ ^0' 



p + 3p + AXp-f + RfR + i3R + eR)fRR 



3 2fn 

+ SR'fRRR^ . (26) 

Combination of Raychaudhuri and Friedmann equations yields 

i?= -2(0 + ^02). (27) 

We assume that there exists an exact power law solution to the modified field 
equations 

a{t) = ao^, (28) 

where m is a positive real number. If < m < 1, then the required power 
law solution is decelerating while for m > 1 it exhibits accelerating behavior. 
For the equation of state p = up, the energy conservation equation leads to 

p{t) = pot-^''<^+'^\ (29) 

Using Eq.f l28|) in Eq. fl27|) . the scalar curvature becomes 

R = -6m(2m - l)r^ = -r]^r^, (30) 

where 77^ = 6m(2m — 1). We see that sign of R depends on the value of m, 
i?>OifO<m<| and i? < for m > |. Since m = ^ leads to vanishing 
of R, so we exclude this value of m in our discussion. 

Using Eqs. fl2^ and fl5U]) . Friedmann equation can be written in terms of 
Ricci scalar R, f and Jr as 

3m(l+cij) 

fRRR' + ^^RfR + ^^^f-{2m-l)ApJ^) ' =0, (31) 
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where A = 1 + ^(3 — u). This represents second order differential equation 
for f{R) whose general solution is 

3m(l+u;) 

f{R) = Xmco (—) ' + Cii?3(3-"»-^) + C'2i?i(3-"^+V3;^), (32) 



where 



y 4A(2m - l)po _ 2, in .1 

3m2(3a; + 4)(a;+l)-m(9a; + 13) + 2' + ium + i, 

and Ci, C2 are arbitrary constants of integration. Since m > 0, so 5^ > 

for cosmologically viable solutions. X,„c^ is found to be real valued but it 
diverges for 3m^(3a; + 4)(w + 1) — m(9uj + 13) + 2 = 0. i.e., m and u satisfy 
any of the relations u = ^'^"^^ or m = ^•^+6^'^+i"j3++i y^- Since i? < 0, so 
{—R/r]m) > for all R, thus we have real valued solution f{R, T) = f{R)+XT 
showing that the power law solution exists for this model. 

For A = 0, we obtain solution as in f{R) gravity^^\ To check whether 
the f{R,T) gravity reduces to GR, we need to put Ci = C2 = A = 0. For 
m = and po = 3(i+ui)^ > ^^^^ theory reduces to GR. We are interested 

to construct the f{R, T) model of the form aR^ + AT. If we put m = 3^^"^^ , 
then f{R) is given by 

f{R) = aU-RT, (33) 

where 

_ 2^-^"3"-^nA(n(4n - 3(1 + a;))^-"(l + o;)^""^ 
"""^ ~ (n^(6uj + 8) - n(9uj + 13) + 3(uj + 1)) ' 

and hence f{R,T) — anw{—RY + AT. This model represents the exact 
Friedmann-like power law solution a oc i3(i+") and the limit n — )> 1 with 
A = leads to GR. For n — 1, our solutions represents ACDM model of the 
form/(i?,T) = i? + AT. 

4.1.1 Phantom Phase Power Law Solution 

We construct the phantom phase power law solution which lead to big rip sin- 
gularity. For this case, the scale factor and Hubble parameter are expressed 
as 

777 

a{t) = ao{ts-tr'^, Hit) 



ts-t 
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The scale factor diverges within finite time {t — )■ ts) leading to big rip sin- 
gularity for m ^ 1^^^ The results for this case can be recovered by just 
replacing m by —m in the previous section. Hence, the phantom phase 
power law solution exist for f{R) + XT gravity. 



4.1.2 Constraining f{R, T) = f{R) + XT Gravity 



In section 2, we have found that Raychaudhuri equations with attractive 
behavior of gravitational interaction give rise to SEC and NEC which hold 
for any theory of gravitation. In this form of f{R, T) gravity, one can employ 
an approach similar to that in GR to develop the energy conditions. 
Equations ([6]) and (122!) can be written as 



R^ 



111/ 



rj. _T 

J- /II/ ~^9/j,i/) 



where 



T, 



T 



1 

Jr 
1 

Ir 



1 + X)T^^ ^\^P~ P'>3/iu + (V/,Vi. - g/cuO)fR 
[(1 + 3A)T + AXp + f- RfR - SD/h] . 



Equations (fT4|) and (134|1 lead to following inequalities 



NEC 
SEC 



1 

Jr 
1 



Vr 

+ RH)fRR 



(P +P){l + X) + {R- RH)fRR + R^fRRR 

(p + 3p) +ApX-f + RfR + SR^fRRR + 3iR 
> 0. 



(34) 



> 0, (35) 



(36) 



For A = 0, one can obtain the NEC and SEC in f{R) gravity. Furthermore, 
the more familiar forms of NEC and SEC in GR, i.e., p+p ^ and p+3j9 ^ 0, 
can be recovered if A = and f{R) = R. 

To derive the WEC and DEC, we can extend the GR approach by intro- 
ducing an effective energy-momentum tensor. The above inequalities of the 
SEC and NEC are obtained directly from Raychaudhuri equations, however 
equivalent results can be derived by using the conditions pe// +Peff ^ and 
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Peff + 'ipeff ^ 0. From Eq. 
pressure are given by 

1 

7^ 



the effective energy density and effective 



Peff 



p+^i^p-p) + kf 



Peff 



+ 



1 

R'f. 



X 



(p - -(p - 3p) - -(/ - Rfn) + {R + 2RH)fRR 



RFiR 



(37) 



(38) 



p + \{^p-p) + \u 



Rfn) - 'iHRfRR 



^ 0, 



(39) 



{p-p) + 2A(p + p) + f-RfR-{R + 5RH)fRn - R?fRRR ^ 0. (40) 



The WEC and DEC in f{R) + AT gravity can be obtained by following 
the effective energy-momentum tensor approach. The WEC is obtained by 
satisfying inequality (135]) and the constraint 

1 

7^ 

and DEC is obtained by satisfying inequalities fl35|) . f p9|) and the constraint 
1 

Ir l 

We find that by setting A = 0, the WEC and DEC in f{R) gravity can be 
recovered. Moreover, for f{R) = R and A = 0, the well-known form of weak 
and dominant energy conditions in GR can be reproduced. 

The above energy conditions can be used to put constraints on a given 
f{R) model in the context of f{R,T) gravity. We assume that > to 
keep the effective gravitational constant positive. Using the relations (!2T]) . 
the energy conditions for f{R) + AT gravity in terms of present day values 
of if, q, j and s are given by 



NEC 



WEC 



X 



(1 + A)(po + Po) - GH\so - Jo + (go + l)(go + S))foRR + H^[6 
HoUo - go - 2)ffoRRR ^ 0, 



A 



1 



Po + 2 (3Po - Po) + + 3i^o (1 - Qo)foR + 18H^[m - Qo - 2] 

/ORR ^ 0, Peff + Peff >0, 

SEC : (po + 3po) + 4poA - /o - 6H^{1 - qo)foR + 3[6iJo'(jo - go - 2)]^ 



X 



X 



DEC 



foRRR - l8H^[so + Jo + go + 7go + A]foRR ^ 0, 
(Po - Po) + 2A(po - Po) + fo + ei^o (1 - go) /or 



Peff + Peff > 0, 

m'oUo - go 



2)?foRRR - GH'lso + 5jo + (go - l)(go + moRR > 0, 

Pe//+Pe//^0, Pe//^0. 
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In order to present the concrete application of the above energy conditions 
in f{R,T) gravity, we employ the exact power law solution of f{R) + XT 
gravity. We consider the present day values of deceleration, jerk and snap 
parameters as go = -0.81 ± 0.14, jo = 2.16l|J:?^ and sq = -0.22l°;?^28)_ ^Ve 
shall discuss the WEC requirement to illustrate how the above conditions 
can be used to place constraints on f{R,T) gravity. We note that all the 
above conditions depend on the present value of pressure Pq, so for simplicity 
we assume p = 0. 

Now, we take the power law solution as an objective model which is given 

by 

/(i?,T)=«„(-i?)" + Ar, (41) 

where n is an integer and a„ = '^(g„2'!^3"^^^_3)^"'' — • The constraints to 

fulfill the WEC, i.e., peff ^ 0, pe// +Peff ^ 0, are respectively obtained as 

(2 + 3A)po + a„[6ifo'(l " QoTlBiin' -n)-n + l]^0, (42) 
(1 + A)po + a„n(n - l)6H^[QH^,il - go)]""'[-(% - Jo + (go + 1) 
X (go + 8))-52(n-2)] ^0, (43) 

where Bi = (jo - go - 2)/(l - qo f and B2 = (jo - go - 2)V(1 - go)- As the 
standard matter is assumed to satisfy the necessary energy conditions and 
A > 0, so (2 + 3A)po > and (1 + A)po > 0. Hence, the inequality is 
reduced to 

a„(3.3ifo)^"/3n ^ 0, where = Bi{n^ - n) - n + 1. 

It is clear from above expression, the result is trivial for n = 0, 1. We consider 
the following two cases: 

(i) an > 0, the allowed values for n are n = {2, 3, 4, ...}. Now /3„ > in the 
range n = {4, 5, 6, ...} and < for n = 2, 3. 

(ii) a„ < 0, the acceptable values of n are n = {—1,-2,...} and in this 
particular range we have /3„ < 0. Thus, the inequality pe// ^ is satisfied 
forn = {...,-2,-1,4,5,...}. 

Now we check the validity of Eq. pSj) except n = 0, 1 as the result is trivial 
for this choice. The inequality is transformed to the following form 

a„(3.3/fo)^"~Vn ^ 0, where p„ = (n^ - n)(2.054 - 0.52n). 

The results of the above inequality can be interpreted as: 

(i) p„ > 0, if n = {2, 3, —1, —2, ...} and for p„ < 0, the acceptable values of 
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n are n = {4, 5, 6, ...}. 

(ii) > 0, with acceptable range n = {2,3,4,...} and a„, < 0, when n = 
{— 1, —2, —3, ...}. Hence, the condition Peff+Peff ^ is satisfied for n = 2,3. 



4.2 f{R,T) = R + 2f{T) 

Now, we construct the power law solutions for R + 2/(T) gravity, where 
/(T) is an arbitrary function of the trace of energy-momentum tensor. The 
effective Einstein field equations are given by Eq.([6]) with 

T^/J = (1 + 2/t)T^. + (2p/t + f)g,., 

where /t is the derivative of / with respect to the trace of energy-momentum 
tensor T. The Friedmann equation and the trace equation can be obtained 

as 

&^ = 3[p + 2ip + p)fT + f], R=-{p-3p)-2{p + p)fT-Af. (44) 

It follows that the field equations can be represented as the Raychaudhuri 
equation 

® + \®^ = -\ [(p + 3p) + 2(p + p)fT - 2/] . (45) 



Combining Eqs. fl44|) and fl45|) . we can get the Ricci scalar R given in 
Eq. fl271) . Using Eq.f l29p . the Friedmann equation can be written in terms of 
trace of the energy-momentum tensor T, /(T) and its derivative with respect 
to T as 

Tf , T , (l-3u;)/ K{1 - 3uj)T^^ 

— 2 

where K = 3m^(po(l— 3a;)) 3™Tt+^. This is the first order differential equation 
in /(T) whose solution is 

/(T) = ^— + L^T^^ + CiT=^, (47) 

iij — o 

-2 

where L^lj = ^"^ ^^"^4+3^(1- 3^;)^ ^^"^ ^1 arbitrary constant of integra- 
tion, Lmuj is finite and real valued unless 4 + 3m(l — 3a;) = 0. In general, the 
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function f(T) is real valued if m and co do not satisfy the relation m = 3(]^_3^-) 
and if w > 3. Therefore, the power law solutions exist for R + 2f{T) gravity. 

For m = 0, we have a = Oq, so that H = R = 0, it represents the Einstein 
static universe. The solution for this case is 

/(T) = ^— + CiTTO^. (48) 

(jJ — o 

The standard Einstein gravity can be recovered for the choice Ci = 0, m = 
g^^:^ and po = 3(1^^- order to develop a more general form of function 
/(T), we put m = g^^^, so that 

/(T) = + a^^T-, (49) 

where anuj = 4(i+t^)+2ra(i-3i^) • ensure that this theory 

reduces to GR for n = 1. It is remarked that phantom power law solutions 
exist for this form of f{R,T) gravity, which can be obtained in a similar 
fashion as in section 4.1.1. 

4.2.1 Constraining f{R, T) =R + 2/(T) Gravity 

The effective energy density pe// and effective pressure Pe// for this particular 
f{R,T) gravity are defined as 

Peff = P + 2(p + p)fT + f, Peff =P- f- (50) 

Using Eq.( l50|) in energy conditions ( fT7|) -(l20l). the following form is obtained^®^ 



(p + p)[1 + 2/t] ^0, 

p + 2{p + p)fT + f ^0, Peff + Peff >0, 

p + 3p + 2(p + p)fT -2/^0, Peff +Peff> 0, 

p-p + 2{p + p)fT + 2f ^0, Peff^O, Peff + Peff ^0. 



NEC 
WEC 
SEC 
DEC 

To check how these conditions place bounds on power law solution fHOl) in 
R + XT gravity we put p = 0, so that T = p. Hence, the function /(p) is of 
the form 

/(p) = -S + «nP", (51) 
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where a„ = ^^^^ ^" • "^^^ constraints to accomplish the above energy 

conditions are obtained as follows: 

NEC : ^ + 2na„p" ^ 0, 

WEC : - + 5na„p" ^ 0, 
3 

SEC : ^ + 2(2n + l)a„p" ^ 0, 

DEC : p+{7n + 2)a„p" ^ 0. 

Above conditions are trivially satisfied for n = 0,1. The quantities na„, 
(2n+l)a„ and (7n+2)a„ are negative when n = {—3, —4, —5, ...} and positive 
for n = {—1, 2, 3, ...}. Since p is assumed to be positive, so it is obvious that 
these conditions are satisfied within the range of n = { — 1, 2, 3, ...}. 

5 Stability of Power Law Solutions 

In this section, we are interested to study the stability of power law solutions 
against linear perturbations in f{R,T) gravity. First, we assume a general 
solution H(t) = Hhit) for the cosmological background of FRW universe 
that satisfies Eqs.( 12^ and (jH]). The matter fluid is assumed to be dust and 
evolution of the matter energy density can be expressed in terms of Hh{t) as 

Ph(t)=Poe-^^^'^(*)'^*, (52) 

where po is an integration constant. Since the matter perturbations also con- 
tribute to the stability, so we introduce perturbations in Hubble parameter 
and energy density to study the perturbation around the arbitrary solution 
Hhit) as follows^o) 

H{t) = Hh{t)il + 6{t)), p{t) = Phil + SUt))- (53) 

In the following, we develop perturbation equations for two specific cases 
fiR, T) = f\R) + AT and /(/?, T) = R + 2/(T). 

5.1 f{R,T) = f{R) + XT 

To study the linear perturbations, we expand function f{R) in powers of Rh 
evaluated at H(t) = Hh(t) as 

f{R) = f^ + f^^{R-R,) + 0', (54) 
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where function f{R) and its derivative are evaluated at Rh- The term O"^ 
includes all the terms proportional to the square or higher powers of R. The 
Ricci scalar R at H{t) = Hh{t) is given by 

RH = -6{Hh + 2Hl). (55) 

By introducing the expressions fl53l) and fl5^ in the FRW equation fl23l) . the 
equation for the perturbation 6{t) becomes 



^HhRhf^n 



6{t) + c{t)6{t) = (56) 



where 



c{t) = [In {H-'Rirnn)] + H,[2d[ln{f'^)])-' - 1]. 



The conservation equation (125!) implies the second perturbation equation as 

Lit) + 3Hh{t)5{t) = 0. (57) 

We can eliminate 6{t) from Eqs. (!56|) and (l57Il and arrive at the following 
second-order perturbation equation 

Lit) + c^{t)Lit) + TTTT^^^- = 0' (58) 

where 

ci(t) = I [/n (i/,-^/2Ui^«)] + H,[2{j^[ln{f^^)])-' - 1]. 

Here, we consider the f{R,T) model proposed in section 4.1 for the dust 
case which is defined as f{R) = an{—R)"' + XT. We evaluate f{R) and its 
derivatives at H{t) = Hhit) and hence the perturbation 5m{t) is given by 

5^{t)=C+t^^ +C^t^-, (59) 

where C± are arbitrary constants and 



_ - ISn + 13 ^ v^n2(8n2 - 15n + 3)2 + 18(8n3 - 2ln^ + 16n - 3)po 
6(r7, — 1) 6?7,(ra — 1) 
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In order to study the stability of perturbation given by Eq.( 1S^ . one needs 
to check the signs of exponents The exponents are found to be negative 
provided that n ^ —2, otherwise fi± would be positive and the perturbation 
is unstable. The perturbation 6{t) is found to be 

m = + (60) 

where = /i± — 1. It can be seen that exponent is negative for n ^ —2 
and is always negative. Hence, as the time evolves the condition n ^ —2 
ensures the decay of perturbations 6{t) and Sm{t) which implies the stability 
of power law solution for this f{R, T) gravity. 



5.2 f{R,T) = R + 2f{T) 

We explore the behavior of perturbations f l53|) for this f{R,T) model and 
expand the function in powers of Th{= Ph) as^^^ 

f{T) = f + f!^{T-n) + 0\ (61) 

where O term includes all the terms proportional to the squares or higher 
powers of T. The function /(T) and its derivatives are evaluated at T = T^. 
Using Eqs.dHS]) and (EI]) in FRW equation (jH]), it follows that 

in + STJ^ + 2T^f^T)6m{t) = QHl6{t). (62) 

Combining Eqs.( l57|) and ( l62l) . the first order matter perturbation equation 
is 

^Ut) + 7^iTh + mfr + 2Tlf!^T)5ra{t) = 0. (63) 

which leads to 

5m{t) = C4exp 1^ y Crdt^ , Ct = ^(1 + 3/^ + 2T,/^^). (64) 
The behavior of perturbation 5{t) can be seen from the relation 
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We explore the stability of power law model (proposed in section 4.2) of 
the form 

/(T) = aiT + asT", (66) 

where oi and 02 are parameters. One can evaluate the expression Ct and 
integral ^ / Crdt for the model flMl) as 

3 

Po(3ai + l)t 



3 



5n n 



1) 



Po(3ai + 1)^2"+! + a2p^n(2n + 1)^2"'+! 

1) 



n + 1 



-2(n2-l) 



,(67) 
.(68) 



As the time evolves, we need to set the conditions for decay of perturbations. 
It is obvious that expression (1671) and ( 168|) decay as time increases for the 
choice n > 1 which results in decay of 5{t) and 6m{t)- Hence, for large 
values of t perturbation decays, this corresponds to the stability of power law 
solutions for R + 2/(T) gravity. We find that the conditions developed for 
stability are compatible with some constraints to fulfil the energy conditions. 
Hence, we may remark that power law solutions are acceptable regarding to 
the stability, energy conditions and late time acceleration of the universe. 



6 Conclusions 

The issue of accelerated expansion of the universe can be explained by taking 
into account the modified theories of gravity such as f{R,T) gravity. The 
f{R,T) gravity provides an alternative way to explain the current cosmic 
acceleration with no need of introducing either the existence of extra spatial 
dimension or an exotic component of DE. In this modified gravity, cosmic 
acceleration may result not only due to geometrical contribution to the total 
cosmic energy density but it also depends on matter contents. This theory 
depends upon matter source term, so each choice of matter Lagrangian C^n 
would generate a specific set of field equations. The various forms of La- 
grangian in this gravity give rise to question how to constrain the f{R,T) 
gravity theories on physical grounds. We have made an attempt to address 
this issue and classify the particular f{R,T) models. In this respect, we 
have developed the energy conditions on general as well as particular forms 
of this gravity. The energy conditions in modified theories of gravity have 
a well defined physical motivation, i.e., Raychaudhuri's equation along with 
attractive nature of gravity. 
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The energy conditions and exact power law solutions are studied in the 
framework of f{R, T) gravity. We have derived the energy conditions directly 
from effective energy-momentum tensor approach under the transformation 
p— j-Pe// and p^Peff- The general inequalities imposed by these conditions 
are presented in terms of deceleration (g), jerk (j) and snap (s) parame- 
ters. In order to get some insights on the application of these conditions, we 
consider two particular forms of f{R,T) gravity, i.e., f{R,T) = f{R) + XT 
and f{R,T) = R + 2/(T). In standard paradigm, the expansion history of 
the universe underwent a power law decelerating phase followed by late time 
acceleration. Therefore, power law solutions are important in cosmology to 
represent the matter dominated phase that later connects to an accelerating 
phase. We have shown that exact power law solutions exist for a special class 
of f{R,T) models. These solutions mimic the ACDM model as particular 
case. We have obtained the necessary constraints to fulfil the energy condi- 
tions for this particular class of f{R,T) gravity. We summarize the results 
of these two models as follows: 

• f{R) + AT 

It is shown that exact power law solution exists for this form of f{R, T) 
gravity given in Eq.f l32p . In the limit of A = 0, the corresponding result can 
be recovered in f{R) gravity. To ensure that this theory reduces to GR, we 
need to set C7i = C2 = A = 0, with m = -^(^^^ and po = 3(1^^- We have 
constructed the general form of f{R) + AT model which corresponds to R"' 
gravity. For this particular model of f{R,T) gravity, the NEC and SEC are 
derived from the Raychaudhuri equation together with the condition that 
gravity is attractive. It is shown that these conditions differ from those de- 
rived in the context of GR and f{R) gravity^^-*. The general expression of 
weak and dominant energy conditions are obtained by introducing the effec- 
tive energy-momentum tensor in the context of GR. We have examined the 
WEC bounds on i?)" + AT model in terms of present day observational 
values Hq, qq, Jq and Sq. 

• R + ViT) 

The model f{R,T) = R + 2f{T) corresponds to gravitational Lagrangian 
with time dependent cosmological constant being function of trace of the 
energy-momentum tensor^^). This model appears to be interesting and has 
widely been studied in literature^"^^-'. A general form of /(T) model (H9|) 
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is obtained which corresponds to GR in the hmit n — 1. Wc have apphed 
the energy conditions to set the possible constraints on this f{R,T) modeL 
It is found that energy conditions are globaUy satisfied within the range 
of n = { — 1, 2, 3, ...}. It is worth mentioning here that results of power law 
solutions and energy conditions obtained in this paper are quite general which 
correspond to GR and f{R) gravity. 

We have also analyzed the stability of power law solutions under linear 
homogeneous perturbations in the FRW background for f{R,T) gravity. In 
particular, perturbations for energy density and Hubble parameter are intro- 
duced which produce linearized perturbed field equations. It is shown that 
stabihty/instabihty can be studied for particular f{R, T) models under some 
restrictions. The stability conditions are found to be compatible with energy 
conditions bounds to some extent. Hence, power law solutions in f{R,T) 
gravity can be considered as viable models to explain the cosmic history of 
the universe. It is also interesting to note that for these f{R,T) models, 
the gravitational coupling becomes an effective and time dependent coupling 
which modifies the gravitational interaction between matter and curvature. 
Our analysis shows that stable power law solutions are contained in class 
of f{R,T) models, at least considering a given background evolution of the 
universe. 
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